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Abstract 

We compute the differential and the total cross sections for pp scattering using 
the QCD pomeron model proposed by Landshoff and Nachtmann. This model is 
quite dependent on the experimental electromagnetic form factor, and it is not 
totally clear why this form factor gives good results even at moderate transferred 
momentum. We exchange the eletromagnetic form factor by the asymptotic QCD 
proton form factor determined by Brodsky and Lepage (BL) plus a prescription 
for its low energy behavior dictated by the existence of a dynamically generated 
gluon mass. We fit the data with this QCD inspired form factor and a value for the 
dynamical gluon mass consistent with the ones determined in the literature. Our 
results also provide a new determination of the proton wave function at the origin, 
which appears in the BL form factor. 
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The increase of hadronic total cross sections was theoretically predicted many 
years ago ^ and this prediction has been accurately verified by experiment At 
present the main theoretical approaches to explain this behavior are the Regge pole 
model and the QCD-inspired models. 

In the Regge pole model the increase of the total cross section with the energy 
is attributed to the exchange of a colorless state having the quantum numbers of 
the vacuum: the pomeron |3|] . According to the Regge pole theory the hadronic 
total cross section at high energy (s ^> \t\), where s and t are the Mandelstam 
variables of the process, behaves as a to t oc s^ ' 1 . a(0) is the intercept of the Regge 
trajectory a(t) = «o + a't, where «o = 1 + £o, eo ~ 0.08 and a' = 0.25 GeV~ 2 . 
This behavior leads to a violation of the Froissart bound, although it is believed 
that the exchange of more pomerons eventually will lead to the unitarisation of the 
scattering amplitude. 

Within Regge theory Donnachie and Landshoff [4| predicted the differential elas- 
tic proton-proton scattering to be given by 

^ = const. xFi^Vs)^- 1 ' , (1) 

where Fi (t) is the proton form factor measured in ep elastic scattering. This expres- 
sion fits the data quite well at small t |5j , and is reasonable even at moderate values 
of t. It is clear that this cross section gives a non-trivial check for the functional 
form of Fi(t), as well as it not clear at all why the electromagnetic form factor works 
so well for the pomeron. 

In this work we will consider a QCD pomeron model to compute Eq.fll}. This 
model has already been used successfully, therefore we use the data in order to 
discuss its dependence on the proton form factor. The idea is to exchange the 



eletromagnetic form factor by t 



re asymptotic QCD proton form factor determined 



by Brodsky and Lepage (BL) f|. Of course this last one is not expected to be 
suitable at the infrared (IR) region, but it is for this region that we introduce the 
concept of a dynamical gluon mass and one hypothesis about how the form factor 
should behave in the IR, matching it to the BL proton form factor at high energy. 
Our model for the proton form factor is used to obtain the differential and the total 
cross-sections for pp scattering, and it fits the data with a value for the dynamical 
gluon mass consistent with the ones determined in the literature. In this calculation 
we obtain a reasonable value for the proton wave function at the origin that appears 
in the BL form factor, which is used in many phenomenological applications, and 
turns out to be related to the dynamical gluon mass. 

In the QCD framework the pomeron can be understood as the exchange of at least 
two gluons in a color singlet state 0]. A model for the Pomeron has been put forward 
by Landshoff and Nachtmann where it is evidenced the importance of the QCD 
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non-perturbative vacuum [8]. One of the aspects of this non-perturbative physics 
appears as an infrared gluon mass scale which regulates the divergent behavior of 
the Pomeron exchange [9|. 

The Landshoff and Nachtmann (LN) model can explain diffractive scattering 
data quite successfully 0, 10], at the same time that it is a quite simple model. The 
calculation of hadronic cross sections in this model are straightforward, although 
there are two approximations that are usually performed in order to compare the 
model to the data. First, the pomeron exchange dependence on the energy has to be 
introduced in an ad hoc way. This means that we multiply the scattering amplitude 



by a factor of the form Z = (s/wq) 



2\0.08+a't 



1/a' is an energy scale. We know 



that such behavior must come from the exchange of multiple gluon ladder diagrams 
described by the Balitski-Fadin-Kuraev-Lipatov (BFKL) equation ll|. However we 
also know that the exponent of the BFKL is not simply related to the data. Of 
course we expect that in the future the LN approach can be made compatible with 
the BFKL one. In the meantime we just follow as described above and introduce 
the factor Z by hand. Secondly, the cross sections depend on the square of hadronic 
wave functions, which are ultimately related to their respective form factors. As 
discussed above, in general it is the electromagnetic phenomenological form factor 
that is used in the actual diffractive cross section calculation ;9']. The ideal case 
would be the use of a proton form factor totally derived from QCD, and this is the 
point that will deserve our attention. 

In the LN model the elastic pp differential cross section can be obtained from 



da 
~dt 



\A(s,t)\'< 



Wits 2 ^ 
where the amplitude for elastic proton-proton scattering via two-gluon exchange can 

be written as 



with 



A(s,t) = is8a 2 s [Ti - T 2 



(3) 



,MD[l + k\D(Z-k)\G l ,(<i.0)\- 



(4) 



T 2 = J d 2 kD (| + ^(|- k ) G P (q,k 



2G p (q,0)-G p (q,k 



where G p (q, k) is a convolution of proton wave functions 



(5) 



G p (q,k) = / d 2 pdKip*(n,p)ip(K,p — k — nq). 



(6) 
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To estimate G p (q, k — q/2) it is usually assumed a proton wave function peaked at 
K — 1/3 implying that 



/, 2 - £ 
4 



(7) 



In the many calculations of this model up to now G p (q, 0) was given by the Dirac 
form factor of the proton 

, , „ , n . Am 2 - 2.79t 1 
F,( f ) = G,(,,0) = 4m2 _ f (1 _ i/Q71)2 . (8) 

This means that the strongly interacting pomeron sees the proton in the same way 
as it is seen by a photon. This is not totally surprising, as discussed in the work of 
Landshoff and Nachtmann [8|, since the pomeron formed by gluons that propagate 
in the vacuum up to a certain critical distance (the gluon mass in our case) lead 
naturally at high energies to a pomeron that interacts as an isoscalar photon. 

In this work we would like to differ from the previous calculations exactly by the 
use of the asymptotic form factor expression determined by Lepage and Brodsky 
[l^j ]. which is given by 

F Q cn(t) = G p (q, 0) « , (9) 

where (3 = ll — (2/3)nfi avar is the coefficient of the QCD f3 function. The g 4 behavior 
comes from the propagators of two gluons exchanged between the quarks forming 
the proton and C is determined by the qqq wave function at the origin. Contrary 
to the electromagnetic form factor given by Eq.flSJ) the one of Eq.Q is just the 
asymptotic one, and not normalized to 1 at q 2 = 0. 

We know that the BL proton form factor should be reliable at asymptotic energies, 
i.e. it is just the tail of the actual form factor. How can we perform a full QCD 
calculation of the differential pp elastic scattering without knowing the IR behavior of 
the form factor? To guess something about its IR functional form we must introduce 
in the form factor calculation the concept of a dynamically generated gluon mass, 
which is already going to modify the expression of Eq.fjHJ). 

Eq.(J2J) up to Eq.(0) were computed in Ref. 0, using nonperturbative gluon 
propagators endowed with a dynamical gluon mass, whose existence also imply in 
an IR finite coupling constant 13j. The possibi lity that QCD generates a dynamical 
gluon mass has been put forward by Cornwall |l4j ]. and there are evidences for such 
behavior obtained from solutions of Schwinger-Dyson equations (SDE) 14, ^ and 
from lattice simulations (see Ref . 16| and references therein) . These propagator and 
coupling constant have been used in many phenomeno logical calculations that are 
sensible to their infrared finite behavior [lO,LL7j, and are given by: 
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<hAv_ 
q 2 



(10) 



where the expression for D(q 2 ), that was obtained by Cornwall SIS db fit to the nu- 
merical solution of a gauge invariant gluonic SDE, is equal to 



q 2 + Ml 



bg 2 In 



q 2 + AM 2 



A 2 



where M g (q 2 ) 



is a dynamical gluon mass described by, 



M, 



ln(- 


l 2 +4m g 2 \ 


A 2 ) 


ln 


( 4m 9 2 \ 
V A 2 ) 



-12/11 



(12) 



A(= Aqcd) is the QCD scale parameter. The infrared finite coupling constant 
obtained in the same procedure has the following expression 0] 



®s{q = r 

/31n [(q 2 + 4M 2 (q 2 ))/A 2 
A typical phenomeno logical value for the dynamical gluon mass is 19. I 



m g = 500 ± 200 MeV 



(13) 



(14) 



for A = 300 MeV. Finally, note that these solutions have been obtained in the 
Euclidean space, as well clS £1X6 the momenta in the above equations. 

In practice the SDE equations have to be solved with approximations and there 
are different functional forms proposed for the gluon propagator. In particular, a 
quite simple expression has been proposed in Ref. [3| 



D-\q 2 



Q 2 + M 2 (q 2 ] 



(15) 



and this expression, even neglecting the q 2 dependence in M 2 (q 2 ), is the one that 
we shall use in our calculation. The same happens to the IR behavior of the run- 

and since the 
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Q, 



ning coupling constant. Its behavior at q 2 = is around 1 
calculation is not strongly dependent on the ultraviolet logarithmic behavior of the 
coupling constant, we will just assume a constant value for this quantity. The effect 
of these approximations will be commented ahead. 

Once we admit the idea of a dynamically generated gluon mass we can guess 
what could be expected for the proton form factor. Asymptotically, if we neglect 
the a s dependence in Eq.@, we have 



G p {qMc 



c 



(q 2 + m 2 ) 2 



(16) 
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It could be argued that this approximation is rough due to the fact that we neglected 
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the dependence on a s . Actually this is not so trivial, because Cornwall 
out that the product g 2 D(q 2 ) is independent of the coupling constant. In Ref.[l7( we 
have shown that for gluon masses in the range of Eq. (jl4|) the values of the infrared 
coupling constant vary between 0.5 and 1. We will arbitrarily use the infrared value 
of 0.8. Unfortunately we still have a poor knowledge of the infrared QCD behavior, 
but if the Cornwall's argument is correct it may happens that our approximation is 
much better than it looks like. 

Since the gluon has its propagation limited by the value of the gluon mass {m g ), 
we expect that the IR proton form factor is approximately constant in a radius 
determined by this mass, and falls off vary fast after that. The simplest choice that 
we can make for this IR behavior is 

G p (g,0)|^«exp[-g 2 /2mJ] , (17) 

which is naturally normalized to 1 at the origin. Note that more complicated ex- 
pressions could be used to describe the IR behavior, but they barely would reflect in 
such a simple way a proton that is formed by the exchange of dynamically massive 
gluons. 

The full form factor must be a combination of Eq. (jl6|) and Eq. ffTTj) . and must 
match at one intermediate scale, which we propose as being the dynamical gluon 
mass scale, with the Brodsky and Lepage proton form factor. Therefore our ansatz 
for the proton form factor for the full range of momenta is given by 

G p (q, 0) = exp[-q 2 /2m 2 g }6(2m 2 g - q 2 ) + - + 3 ^ 2)2 % 2 - 2mJ) . (18) 

Eq. (JT8j) is correctly normalized at the origin, match softly the low and high energy 
of the proton form factor, and, if compared to Eq. ([TK]) . gives a prediction for the 
proton wave function at the origin 

C=— £, 19 

e 

whose value will be determined when we compare our cross section calculation to 
the experimental data. 



We compute t 



re differential pp scattering cross section in the LN model, as 



o 



Eq. ipjl . In Fig.(l) 



performed in Ref . j| , with the help of the proton form factor 

the result is compared to the experimental data from et. al. jjjj]. The data are well 
fitted assuming a dynamical gluon mass of m g = 460 MeV, for A = 300 MeV. 
The calculations have some dependence on the gluon mass. We also present in Fig. 
(1) the results obtained using m g = 400 and 550 MeV. We also compute the total 
cross section obtained from Eq.Q (after multiplication by the factor Z). It is clear 
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that the LN model calculation, which turns out to be a function of m g and C, can 
reproduce only the Pomeron term (the one that increases with s) of the fit for the 
total cross sections proposed by Donnachie and Landshoff 0], where it is shown 
that the high energy behavior of pp (pp) scattering is proportional to 21.7s 0808 . 



is propoi 

The high energy data on pp scattering from Ref. |2l| is also shown in Fig. (2) 
and compared with our calculations, using the same values of the dynamical gluon 
mass, which turns out to be, together with the Donnachie and Landshoff fit, just 
another check for our results. 

From Fig. (1) and (2) we see that the data shows a preferred value m g = 460 
MeV, that is compatible with previous determinations of the dynamical gluon mass. 
With this value we obtain 

C = 0.15 GeV 4 . (20) 

It is interesting that in Ref. 0] the data of elastic proton-proton scattering is well 
fitted assuming a dynamical gluon mass of m g = 370 MeV for A = 300 Me^, 
while our recent calculation in Ref. [2^ shows excellent agreement with the data for 
m g = 400 (—100) (+350) Mel^. Our result is in agreement with such calculations 
and provides a new determination of the proton wave function at the origin. 

The value of C is quite compatible with the results obtained from the branching 
ratio tv r%^r^ s 1 whose experimental value is 0.0022. In this calculation we use the 
expression for the branching ratio given by |2 

T(lp-^pp) 001n6 3/ \ \Pc.m.\ < T > 2 



3.210X(s) ' ^' 4 (21) 



r(^ — > hadrons) s y/s s 
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Fig. (2) Proton-proton total cross section. The experimental data {pp) are from Ref. 21 1 
It is also shown the curve given by the Donnachie and Landshoff fit. 



where \P c . m .\/Vs ~ 0.4, s = 9.6 GeV 2 and 

<T>= / [dx\[dy\ - — r \Tt — s 

Jo ViVm Fi(l -Vx) + -xi)j[x 3 (l -y 3 ) + y 3 {l - x 3 )\ x l x 2 x 3 

(22) 

We performed the above calculation using 

/ g 2X -2/3/3 

4>(x) = 4>(x) ass = Cx x x 2 x 3 y°9-^2 I (23) 

Thus we find C ~ 0.13 GeV 4 . This value is consistent with the one of Eq.()20|). and 
shows that our simple approach to the proton form factor is quite sensible. 

In conclusion, we introduced in the proton-proton elastic scattering calculation 
within the LN model the asymptotic behavior of the QCD proton form factor, sup- 
plied by a prescription about its infrared behavior in terms of a dynamically gen- 
erated gluon mass. The high energy behavior of the pp diffractive cross section is 
determined as a function of the factor C, which is related to the qqq wave function 
at the origin, and turns out to be obtained as a function of the gluon mass. We 
obtained C = 0.15 GeV 4 after fitting the elastic pp scattering data with m g = 460 
MeV for A = 300 MeV. Our result provide a desirable connection of the data with 
the perturbative calculation of the proton form factor. Besides that we believe that 
the form factor expression presented in Eq. (jl8|) can be further studied in different 
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processes. Finally, the approximations that we have performed are related to the 
poor knowledge of the gluon propagator and coupling constant infrared behavior, 
and any improvement in the calculation could be only performed after an equivalent 
improvement in the expression of the QCD infrared Green functions. 
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